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Oceanic pressure measurements, even in very deep water, and atmospheric pressure or 
seismic records, from anywhere on Earth, contain noise with dominant periods between 
3 and 10 seconds, that is believed to be excited by ocean surface gravity waves. Most 
of this noise is explained by a nonlinear wave-wave interaction mechanism, and takes 
the form of surface gravity waves, acoustic or seismic waves. Previous theoretical works 
on seismic noise focused on surface (Rayleigh) waves, and did not consider finite depth 
effects on the generating wave kinematics. These finite depth effects are introduced here, 
which requires the consideration of the direct wave-induced pressure at the ocean bottom, 
a contribution previously overlooked in the context of seismic noise. That contribution 
can lead to a considerable reduction of the seismic noise source, which is particularly 
relevant for noise periods larger than 10 s. The theory is applied to acoustic waves in the 
atmosphere, extending previous theories that were limited to vertical propagation only. 
Finally, the noise generation theory is also extended beyond the domain of Rayleigh 
waves, giving the first quantitative expression for sources of seismic body waves. In the 
limit of slow phase speeds in the ocean wave forcing, the known and well- verified gravity 
wave result is obtained, which was previously derived for an incompressible ocean. The 
noise source of acoustic, acoustic-gravity and seismic modes are given by a mode-specific 
amplification of the same wave- induced pressure field near the zero wavenumber. 
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1. Introduction 

Ocean waves generate noise in a wide range of acoustic frequencies f s . The upper end 
of the spectrum, f s > 100 Hz, is dominated by wave breaking and associated bubbles 
( Knudsen et al. 19481 ). whereas the lower frequency part, nominally f s <2 Hz is mostly 
expected to be caused by the nonlinearity of the hydrodynamic equations, on which we 
focus here. The general sound generatio n by fluid flows was described theoretically by 
iLighthil] (1952). Longuet-Higgind (|l950l ) showed how seismic waves can be generated by 
the same process, with noise radiating along the Ear th's crust in the fo rm of Rayleigh 
waves. That theory was extended to random waves by iHasselmannl (1 19631). and later cast 
in the more general framework of wave- wave interactions (|Hasselmannlll966l) . Work on 
comp r essibl e flows has also been extended to the study of tsunamis. In that context, 



Okall (|1988f ) discussed the compressibility effect on gravity modes, which we will call 
here 'acoustic-gravity modes', and the gravity effect in seismic 'pseudo- Rayleigh' waves, 
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that we will refer to as 'Rayleigh' modes. Interest in seismic noise has risen sharply over 
the last few years with the en forcement of the Com prehensive Nuclear Test Ban Treaty, 
and motivated by the work of lShapiro et all (|2005l ) who demonstrated that seismic noise 
correlation could provide a unique monitoring method for the properties of the solid 
Earth. 

Recent numerical models based on the Longuet-Higgins-Hasselmann (LH-H) theory for 
Rayleigh wave gen eration have shown good agreement of m odeled seismic noise spectra 
with observations ( Kedar et al. l2008t 1 Ardhuin et ctZ.I [20lTh . It is still unclear whether 
most of the uncertainties on the modeled noise level can be attributed to errors in the 
seismic sources, associated with a poorly constrained directional distribution of the ocean 
surface wave spectrum, or to e rrors in the seismic propagation. Seismic observations hav e 
also revealed body waves (e.g. Koper et all 2010t Landes et all 2010 : rimers et al. 2012 ). 
for which no complete theory has been proposed to date. Vinnikl ( 19731 ) did propose a 
theory for compressional (P) waves, but he did not consider the important effect of the 
water layer. 

Further, the generation of Love waves, which are surface shear waves polarized in the 
horizontal direction, is not well understood. These Love waves are particularly impor- 
tant for frequencies below 0.02 Hz or above 1 Hz (e .g. Bonnefov-Claudet et al. 20061 : 



Nishida et all 120081 : iKurrle fc Widmer-Schnidrig||2008l) . The low frequency Love waves 



may be excited by the direct action of long surfa ce gravity waves (known as infragravity 



waves), on a sloping bottom (IFukao et aLll2010h. in a way similar to the generation of 



primary microseisms described by Hasselmannl |l963 ). Here we narrow the scope of 
investigation and only consider non-linear wave-wave interactions, which leads to noise 
with frequencies double that of the surface gravity waves. 

Finally, the level of acoustic noise has also been explored as a potential s ource of in- 



form at ion on the poorly known dire ction al spectrum of short grav ity waves (jTvler et al 



Il974f ). iFarrell fc Munkl (j2008l l201oh and iDuennebier et all (|2012f ) showed a large vari- 
ability of the spectral level in the frequency range 0.1-50 Hz that is clearly related to the 
sea s tate. Their interpretation of the data, following previo us studies of underwater noise 
(e.g. lHugheslll976l: iLlovdl l ll98lHKibblewhite fc Ewans1ll985l ). is based on sound generated 
by waves in an unbounded ocean. Although this simplified approach is reasonable for high 
frequency noise, the neglected reflection from the seafloor and subsequent reverberations 
in the bounded ocean may strongly amplify the lower frequency resonant modes. 

Given the renewal of interest in seismic and acoustic noise, we found it appropriate to 
revisit Hasselmann's theory. We thus illustrate, correct, and add a few missing aspects. 
These corrections include important terms for intermediate and shallow water that have 
not been considered befor e for the compressible c onditions, although they were verified in 
the incompressible limit (|Herbers &: Guzalll994l ). The compressible equations of motion 
are used to derive a consistent solution in terms of gravity, acoustic and seismic modes, 
including both surface Rayleigh waves and compressional (P) or shear (S) body waves. 
This solution for body waves has not been given before. 

Building on the seismo-acoustic paradigm proposed by lArrowsmith et all ( 2010t ). the 
basic idea of the present paper is that all modes of motion can be excited by ocean surface 
gravity waves (OSGW) of any frequency. For a given pair of interacting frequencies / and 
/', the frequency of the generated noise is / + /', and the different types of waves are only 
distinguished by their phase velocity, or equivalently by their horizontal wavenumber K 
which is the norm of the vector sum of the wave numbers, K = k + k' of the interacting 
OSWG, as shown in figure [D This t ype of wave- wave interaction is one of the lowest 
order interactions ( Hasselmann 19661 ). 

In the physical space, one such interaction excites waves with horizontal wavelength 



double-frequency noise generation by surface gravity waves waves 3 

(a) 

1 10 horizonta 
. . I , i « ■ Phase 

AG a w R P S «c P+S speed (km/s) 




frequency (Hz): 0.02 0.25 0.4 
Wavelength at transition (km) 

S-P+S : 280 22 14 

S-R : 160 13 8 

R - AG : 75 6 4 



Figure 1. (a) The vertical evanescent or propagating nature of the noise field in the solid and 
liquid layers is defined by the horizontal phase speeds relative to the distinct values of the sound 
speed in the ocean (a w ), and the shear (/3) and compression (a c ) speeds in the crust. From 
slow to fast, there are the acoustic-gravity (AG) domain, the Rayleigh (R) wave domain, and 
two body wave domains (5* only, and P and S together) (b) For any fixed frequency, the four 
domains correspond to 4 concentric regions in the wavenumber plane. For three selected noise 
frequencies generated by OSGW in the infragravity, dominant and high frequency ranges of 
the forcing wave field, the limiting wavelengths between the four domains are indicated, using 
a w = 1.5 km/s, /3 — 3.2 km/s, a c — 5.54 km/s. One example of interaction (black vectors) is 
shown with two gravity wave modes that interact to generate a Rayleigh wave. 

Lh = 2n/K, and different vertical patterns in the atmosphere, ocean and crust, due to 
the very different speeds for compression waves in these three media. Figure [2] shows 
how the same forcing can give almost vertical-propagating waves in the atmosphere, 
waves that propagate almost horizontally in the ocean, and evanescent waves in the 
crust. The analysis of a pair of interacting wave-trains is key to our interpretation of the 
different noise modes. The broad-banded wave spectrum of ocean waves results in the 
superposition of all possible pairs of OSGW waves trains, and thus all possible noise waves 
propagating or evanesc ent, that radiate in all directions. This broad spectrum allowed 
lLonguet-HigginsI (|1950l ) to consider all interactions at once, replacing the wave-induced 
forcing by an equivalent point source exerted on the sea surface. However, that latter 
approach is only valid for noise wavelengths much larger than those of the interacting 
OSWG. Table 1 points directly to the main theoretical results of the present paper, and 
all symbols and notations are listed in tables EHU 



2. Water waves theory and noise sources 

Here we give only a brief derivation of t he solution, which i s a straightforward com- 
pressible extension to the solution given by iHasselmannl (| 1962T) . Water column motions 
are expanded in powers of the sea surface slope with a l inear motion for which com- 
pressible effects may be neglected (|Longuet-Higgin s 1950), and a second order motion 
with pressure P2 and velocity potential <f>2- It is for this second order motion that we are 
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shear and compression waves 



Figure 2. Left, schematic of second-order pressure field, in colours, associated with double-fre- 
quency noise, forced by the interaction of a single pair of directionally opposing monochromatic 
wave trains on the sea surface that form periodic groups. The individual waves in the group 
advance very slowly, at the OSWG phase speed. However, in this particular case, the group as 
a whole travels slightly faster than the sound speed in water, and slower than the crust shear 
and compression velocity. As a result, the crust elastic waves are evanescent while ocean and 
atmospheric waves propagate also in the vertical. In the ocean, the superposition of upward and 
downward waves (dashed lines) gives the vertical mode structure. The number of waves in the 
group was reduced for visibility, and for the same reason the amplitudes of the sea surface and 
bottom elevations have been exaggerated as well as the pressure fluctuation in the atmosphere 
relative to those in the ocean. Right, schematic of ocean waves with a relatively broad spec- 
trum, giving rise to the interaction of all possible pairs of wave trains and noise radiation in all 
directions. 



Table 1. Summary of theoretical expressions for the gravity, acoustic and seismic noise fields. 



Surface forcing field 
Bottom forcing field 
gravity-like modes 
acoustic spectrum 
microbaroms 
seismic source (R) 
seismic source (P) 
seismic source (S) 
seismic spectrum (R) 
seismic spectrum (P) 



infinite water depth 
eq. <f2T29|) 
n.a. 



eq. 
eq. 
eq. 



( 3.11 1 



( 3.10 1 



( 3.19 1 



n.a. 
n.a. 
n.a. 
n.a. 
n.a. 



finite water depth h 
eq. ([2728]) 
eq. dpi I 
not given 
eq. (gm I 
eq. (|3~T9 l 
eq. (11732 1 
eq. (I4~35 l 
eq. (14740 1 
eq. gTJSl 
eq. (14745 1 



extending previous results. Following again Longuet-Higgins (1950) we will see that 



approximately obeys the same linear acoustic wave equation as <pi , so that t he forcing of 
the do uble- frequency noise is only due to boundary conditions. In contrast to lHasselmannl 



(1963), we consider here the general finite depth expression for the surface pressure forc- 



ing P2,surf, and in the boundary condition at the ocean bottom the additional forcing 
Vi bot that accounts for the Bernoulli effect of the near-bed orbital wave motion. The 



theory presented here extends the second order finite depth theory of lHasselmannl (| 1962T) 
to a compressible ocean and elastic seafloor. 

2.1. Equations of motion 

We decompose the water density into a mean value p w and a perturbation p <C p w - Ne- 
glecting stratification effects due to temperature and salinity, the fluctuations in pressure 
p and water density p are related by an equation of state, which involves the speed of 
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sound in water a w ( Lig hthilllll978l eq. 32). 

dp = 2^P 
dt Uw dt ' 



(2.1) 



We assume that the motion is irrotational so that the velocity field is given by the 
gradient of the velocity potential 4>. This assumption, which implies a frictionless interior 
ocean, is well supported by the ob served weak attenuation of swells propagating across 
ocean basins ( Ardhuin et al. 20091 ) and loc al comparisons of th e observed wave orbital 
motion with second order wave theory (e.g. iHerbers et al.lll992l ). Vortic ity effects can be 
impo rtant in the vicinity of the surface and botto m boundary layer s (e.g. Longuet-Higginsl 
1970). or in the presence of sheared currents (e.g. |Peregrinel[l976l ). The former effect has 
little influence on the pressure field, which is our primary interest, and we will not consider 
here the effects of currents. Because of the importance of the apparent gravity acceleration 
g, which defines the vertical axis, we separate horizontal and vertical components using 
vectors and gradient operators in the horizontal plane, e.g. u = V</> = (dcf>/dx,d<j>/dy) 
and w — d(j)/dz. The conservation of mass of sea water is 



dp „ 2 
dt = " PwV < 



•dz 



Equations (|2.2j) - (|2.1j) can be combined to eliminate p, 



dp _ _ 2 
dt ~ Pw w 



dz 2 



(2.2) 



(2.3) 



From eq. (|2.1[) . the water density is only a function of the pressure. The two unknowns 
p and <j) are also related by the moment um conserva tion equation, with can be cast in 
the form of Bernoulli's equation (see e.g. lLamblll932L section 20), 



dj)_ 
~dt 



chp 
dz 



p_ 

Pw 



- gz + c(t), 



(2.4) 



with C(t) a time- varying but spatially uniform function. 

The boundary conditions at the surface z = £ are given by the continuity of pressure 
and vertical velocity 







P = Pa 

<h 

= Vc 



VC- 



dt 



at 



c, 



(2.5) 
(2.6) 



with the atmospheric pressure p a . This expression is translated to the mean sea level 
z = with a Taylor expansion for <f> 



dC 
dt 



-V<f>- VC + C 



d 2 ;. 



at 



(2.7) 



Following lLonguet-Higgins (|l950h . we shall now expand the solution in powers of the 
surface slope, with the sea surface elevation d, associ ated with linear w aves, and a non- 
linear correction £ 2 such that \($\ <C |Ci| (see also e.g. lHasselmannlll962[ ). 



C = Ci + C2 , 4> = 4>i + 4>2- 



(2.80,6) 
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2.2. Linear solution 



We consider the case of a constant depth h. Compressibil ity effects in the linear solution 
are negligible for our purpose (see Longuet-Higginsl l"l 950L eq. 123), so that we may use 



Ci 



k.s 



(k«x— sat) 



^ _jfg C ° Sh( „t + f k) ZjJ^-^ (2.90,6) 



k.s 



cosh(fc/i) 



where k is the norm of the horizontal wavenumber vector k, s is a sign index equal to 
-1 or 1, so that s — 1 corresponds to waves propagating in the direction of the vector k, 
and s = — 1 corresponds to the opposit e direction. The radian frequency a is given by 
the dispersion relation for linear waves (|de LaplaceJll776l ) , 



ij gkta,nh(kh), 



giving the group speed 



da a 

9= dk = k 



kh 



sinh(2fc/i) 



(2.10) 



(2.11) 



2.3. The negligible near-surface forcing 
Eliminating p between (|2.3I) and (12.41) , we obtain the acoustic wave equation 

2N 



V 2 + 



dz 2 



d 
If 



d( t> 1 / IT, 



dz 



(2.12) 



where the gravity term has been removed by our approximation of a constant mean 
density. That term is usually neglected eventually in the solution (IStonelevlll926l). The 
expression for the C(t) term in eq. (|2.4I) . is given by eq. (28) in ILonguet- Higgins (1950) 
and it is actually zero provided that there are no standing waves, i.e. we do not consider 
the case of waves trains of exactly equal frequency and opposite direction. For broad 
wave spectra we may indeed neglect these contributions since the measure of such pairs 
of wave components is zero while the measur e of nearly oppos i te way es is finite. 

Our eq. (|2.12[) corresponds to eq. (130) in iLonguet-Higgins (1950), who showed that 
the non-linear terms yield a contribution (his F terms) that is of the order of g/(2a 2 u k) 
times the other terms. Even for very long waves with a wavelength of 50 km, this factor 
is 0.02, so that we shall neglect this source of wave forcing. This leads to a linear wave 
equation for the second order velocity potential 



d 2 c 



dH 



V 2 



c)z : 



for 



h < z € 



(2.13) 



2.4. Surface forcing 

Because the compressibility only affects the mass conservation equation, it does not 
modif y the kinematic and dynamic boundary conditions. These are given bv lHasselmann 
(|l962l) . The unknown £ is eliminated from the linear terms by adding <9 (|2.4j) /<9t evaluated 
with (|2.5I) , and g x (|2.7I) . This combination of kinetic and dynamic boundary condition 
gives an equation for the velocity potential to second order in the wave slope, valid at 
z = 0. Keeping the lowest order non-linear terms we have, 



dt 2 





oz 



V<p! ■ VCi + Ci 



d 2 z 



ld_ 

2di 



50i 

dz 



d 2 cj>t 
dtdz 
(2.14) 
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We now g i ve th e explicit form of the right hand side using eqs. (|2.9 
Hasselmann ( 19631 ). we may rewrite (|2.14j) as 



dt 2 



9JT 
oz 



J2,surf 



Pw 



dt 



As noted by 



(2.15) 



Namely, our problem is equivalent to the effect of a pressure field P2,surf applied at z = 0. 
Using the linear solution 



P2,surf = Pw ^ Dz ( k ' s ' k '' s ') Z t,k Z iM' 

k,s,k' ,s r 



,ie(k,k', 



(2.16) 



with the phase function of interacting wave trains defined by 

9(k, k', s, a') = [(k + k') • x - {so- + s'a') t] (2.17) 

and coupling coefficient D z (k, s, k', s') given by eq. (|A 1[) . 

The shape of this pressure pattern is well understood when the full sum in eq. (|2.16p 
is reduced to only two interacting deep water wave trains with amplitudes a and a' and 
slightly different frequencies a and a', traveling in opposite directions, as shown in figure 
[3] For » 1 we may use D z (k, 1, -k, 1) = -2a 2 and D 2 (k, l,k, 1) = 0. Defining 
K = k — k! , the equivalent surface pressure is given by the sum interaction, 



P2,t 



-2p w a 2 aa {cos [Kx + 2at}} . 



(2.18) 



The corresponding difference interac tion yields a short wave, with wavenumber 2k, that 
will not be considered here (see e.g. iHasselmannl ll 963) . We note that the wave forcing 
P2,surf is out of phase with —u 2 /2, the part of the Bernoulli head which comes from the 
horizontal velocities. Indeed, as the wave trains propagate in opposite directions their 
velocity partially cancel where the surface elevations add up, thus causing a pattern of 
higher pressure under the groups of high waves and lower pressure under the lulls, that 
is opp osite to the familiar Bernoulli eff ect of set-down under groups of uni-directional 
waves (jLonguet-Higgins fc Stewardll962l) . As a result the double- frequency perturbations 
of —u 2 /2 have a sign opposite to P2,surf- That effect will be very important for waves in 
shallow water. In deep water, the contribution of —w 2 /2 must also be considered. Using 
a ~ a' ', we have 



/2- 



(2.19) 
(2.20) 
(2.21) 



u = aa cos(kx — at) — a 1 a 1 cos(kx + at) 
w = aa sin(fca; — at) — a' a' sin(fca; + at) 
-(u 2 + w 2 )/2 = aa'a 2 cos [Kx + 2at] 

Eq. (|2.18|) generalizes the result g iven by Loneuet-Higginsl ( 1950l ) for equal wave pe- 
riods. The standing wave studied bv lLonguet-Higgi ~ ( 1950| ) is thus obtained, somehow 
paradoxically, as the limit of wave groups that travel at an infinite speed, but that are 
infinitely long. 

For wave directions nearly opposite, instead of exactly opposite, this first term can 
propagate in any horizontal direction, given by the direction of K = k + k'. 

2.5. Noise spectrum and finite depth effects 
Both LH50 and H63 used the deep-water approximation 

2" 



P2,: 



+ 



valid for kh ^> 1, 



(2.22) 
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T 2 = 6 s, C 3 = 2000 m/s T' = 11.88 s 




Figure 3. Generation of supersonic wave groups (in blue) by the superposition of two opposing 
deep water monochromatic wave trains of nearly equal periods T and T' . The curves show the 
surface elevation of the individual wave trains (black and red) or their combination (blue). The 
group with length 12 km propagates in the same direction as the wave train with the shortest 
wavelength. If the wave trains are not exactly directionally opposing, the group propagates in 
the direction of the vector K = k + k'. The red and black dots are attached to the wave crest 
of each train, and move 100 times slower than the blue group. 



instead of the more complex but more general form given by eq. (|2.16[) . We will illustrate 
the important differences between these two expressions for kh < 1, by considering the 
forcing of very long components. Indeed, sound waves in the ocean have velocities in 
excess of 1.4 km s _1 . As a result, regardless of the vertical wavenumber I, the horizontal 
wavenumber vector is relatively small K = yflu)* / — I 2 ) < uj s /a w . Thus the acous- 
tically noisy wave interactions verify K -C k, which gives k' ~ — k, and / ~ /' with 
f s ~ 2/. We may thus focus on the estimation of the spectrum of p2, surf at K ~ 0. 

For this purpose we introduce the spectral density of this surface pressure in the three 
spectral dimensions (K x , K y , f s ). Using the Fourier amplitude P2,surf(K, / s ) of the forcing 
pressure P2,surf, with wavenumber vector K and frequency f s , 



F, 



p2,surf 



(K,/ s ) 



2 lim 



|P2,surf(K, f s ) 

AK x dK v df s 



(2.23) 



The factor 2 in the expression makes this a single-sided spectrum, with non-zero values 
only for f s ^ 0. This spectral density is in a three-dimension al spectral space, with S.I. 
units of Pa 2 m 2 s; it is denoted F p3 u in lArdhuin et al. ( 2011 ). 

Using (|2.16l) the surface pressure spectrum can be expressed in terms of quadratic 
products of the (linear) sea surface elevation spectrum 



E(k x , k y ) = 2 lim 



o dk. r dk„. 



with a coupling coefficient from (|A ip that simplifies for K 



to 



D z (k,l,-k,l) = -2 ( x 2 



1 



1 



4sinlr (fc/i) 



(2.24) 



(2.25) 



To transform the spectra to frequency-direction space we use the Jacobian transfer- 
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mation 

E(f,e) = ^E(k x ,k y ). 

We now introduce the directional distribution M such that E(f, 
we define the directional integral 



/(/)= / M(f,0)M(f,6 + n)M 
Jo 

With these notations we finally obtain 

/C 2 dk dk 
E(f,6)E(f,0 + ir) 



(2.26) 

= E(f)M(f,0), and 
(2.27) 



p 2 w g 2 f s E 2 (/)/(/) tanh 2 (^) 



2kh 



sinh(2fc/i) 



1 



i 2 



4sinli (fc/i) 



(2.28) 



In deep w ater (kh ^> 1) , the equivalent surface pr essure is P2,surf — Pw(u 2 + w 2 ), and its 
spectrum (Hasselmann 1963; Ardhuin et all [2011 ) 



F p2 , sur{ (K ~ 0,/ s ) = P 2 w g 2 f s E 2 {f)I{f). 



(2.29) 



As shown in figure [4j this spectral density of surface pressure is, in finite depth, up 
to four times s malle r tha n the deep water approximation obtained from eq. (12.22[) and 
used bv lWebbl (|2007l ) and lTanimotol (|2007t ). This is because the surface pressure actually 
combines two terms. One is p w (u 2 + u;f)/2, from the momentum conservation equation 
- which here takes the form of the Bernoulli equation (12.41) - and the other, equal in 
magnitude in deep water, comes from the non-linearity of the surface boundary condition 
(|2.7p . As kh goes to zero, both the latter term and the vertical velocity contribution to the 
Bernoulli pressure are small compared with p w u 2 /2, thus reducing the surface pressure 
forcing by a factor four. 

This effect is not important for the dominant microseismic peak, generated by waves 
of period about 10 s, but it may be important for the much longer waves, known as hum, 
driven by long surface gravity waves. However, in that case, one should also consider the 
direct action of the wave-induced pressure on the bottom. 



2.6. Additional bottom forcing 

To evaluate noise in the water column, both the second order pressure forcing applied 
at the surface and bottom have to be taken into account. This includes the well known 
Bernoulli effect of a pressure drop in response to an increase in velocity. At the bottom, 
this wave-induced pressure is —p w u\/2 which is exactly out of phase with i?2,surf- In the 
limit kh — > 0, it cancels the source of noise that would have resulted from the surface 
forcing alone. In physical terms, for kh — 0, we have the same momentum balance at 
the sea surface and bottom, and thus the same pressure perturbations, which are zero 
because the pressure just below the surface has to match the atmospheric pressure. 
For any value of kh, the coupling coefficient given by eq. (|2.25l) d iffers from the fu ll 



second order coefficient for the bottom pressure (e.g. eq. 4 in iHerbers fc Guzal 119911 ). 
which also involves the Bernoulli head (the bracket in eq. I2.4[) . However, that extra 
term is also relevant to the generation of seismic noise due to the bottom boundary 
condition that couples the solid crust to the water column. Indeed, the second-order 
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pressure perturbation at the bottom writes, 

Pa(-fr) = ~P™~~^f + P2 : bot, (2.30) 

where the Bernoulli head contribution to the pressure can be expressed from the first 
order wave amplitudes, 

P2,bot = p w D p» ( k < s ' k '' s '< z = -V Z :.k^i k"' ! " k k ' • (2-31) 

k,s,k' ,s' 

with a coupling coefficient D p t given by eq. (|A2|) . 

We may interpret the bottom pressure (|2.30l) as the sum of the surface forcing P2,surf 
transmitted to the bottom by fa, and a direct effect of the Bernoulli head at the bottom 
which is an additional forcing p2,bot that partly cancels P2,surf • 

We shall see in the next section that the forcing term for seismic noise is P2.surf + 
cos(Z/i)p2,bot) with I ^ K <C k the vertical wavenumber in the water. For shallow water 
gravity waves, kh <C 1 and thus cos(lh) ~ 1 so that the effective forcing term becomes 
P2,surf +P2.bot, which equals the bottom pressure in the incompressible limit. The shallow 
water asymptote of the spectrum of this total forcing term is very different from the 
surface pressure only. Compared to eq. (|2.28[) . the [l + 0.25/ sinh 2 (kh)] factor is now 
replaced by 1. For kh -C 1, this is a factor (kh) 4 /16 smaller, as shown in figure |31 
This asymptote is relevant for the hum, the noise with periods large r than 30 s, which 
is believed to be driven by long (infragravity) surface gravity waves (IWebbll2007l) . The 



source of this hum is attenuated by several orders of magnitude on the continental shelves 
and not amplified according to the deep water approximation given by eq. (I2.22p . 



3. From surface pressure spectrum to noise spectra 

3.1. Dispersion relation and modes 

The problem of noise generation by waves has been reduced to that of noise generation 
by an equivalent surface pressure field P2,surf (x, y, t), with the possible addition of a 
bottom pressure field P2,bot(x, y, t) for finite depths. From a statistical point of view, this 
equivalent surface pressure is fully represented by its spectrum, -F P 2,surf (K, f s ) where K = 
k+k' is the sum of the interacting wavenumbers and f s = f±f is the sum or difference of 
the interacting frequencies. Subjected to this surface forcing, our linearized wave equation 
(|2.13[) will have linear solutions. In particular, any propagating or evanescent solution 
will take the following form 

fa oc exp[i(K x x + K y y + Iz - u) a t)], (3.1) 

which, substituted in the wave equation, gives the dispersion relation 

{-uj 2 s +al[K 2 + l 2 ]}=0, (3.2) 

where lo s = 2irf s . Both lu s and K = |K| are imposed by the forcing, so that the magnitude 
of the complex vertical wavenumber I is given by 

which yields 

fa = (Cc a{z+h ^ + Dc- tl< - z+h ^ e ie , 6 = (K x x + K y y - u.t) (3Aa,b) 
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Figure 4. Compared to deep water, the source of seismic noise power in finite water depth is 
amplified by a factor tanh 2 (fc/i) [1 + 2fc/i/sinh(2fc/i)], when accounting for both the surface and 
bottom forcing. This is very different from the approximation that considers the surface forcing 
only (leading to eq. I2.28|> . even more so when using its deep water approximation (eq. I2.22|l 
which gives an additional factor 4 difference for kh <gil. 



where C and D are the bottom amplitudes of the upward and downward propagating 
waves, determined by the surface and bottom boundary conditions. 

There are two classes of solutions. Those for which uj s /K ^ a w and thus I is imaginary: 
these are 'acoustic-gravity ' modes with an amplitude that decays exponentially from the 
surface. For the shortest components, we have uj s /K <C a w and thus I — ifc, corresponding 
to the incompressible limit in which only gravity is important. The other class of solutions, 
for which oj s /K ^ a w and I is real, are the acoustic modes that propagate along both 
vertical and horizontal dimensions, as illustrated in figure [5] 

The equations of elasticity in the crust yield the same wave equation, one for com- 
pressional and another for shear motions, with a w replace d by the crusts compre ssional 
wave speed a c and shear wave speed /3 S , respectively (e.g. lAki fc Richards! 120021) . Since 
we have a w < j3 s < a c , we can distinguish four different regimes depending on the value 
of the horizontal phase speed ui s /K, relative to these three velocities (fig. [TJ. 

In the family of gravity-acoustic modes, there is no resonance in the forcing, namely 
there are no free waves with K = k + k' a nd f s = f ± f and thus the motion is locked 
with the forcing wave groups, as verified bv lrierbers fc Guzal ( 19941 ). The acoustic modes 
propagate obliquely to the vertical away from the surface, and interact with the ocean 
bottom (figure [5|) . The seismic Rayleigh waves modes propagate horizontally, combining 
an acoustic- like motion in the water and evanescent elastic waves in the crust, illustrated 
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L= 2 % [K 




Figure 5. Schematic of wave groups and forced acoustic and seismic wave motion. For a given 
wave group period T, the horizontal wavelength L — 2n/K can be larger than the acoustic 
wavelength in the water due to the oblique incidence of the sound waves. The superposition of two 
obliquely propagating sound waves (arrows) forms a mode pattern that propagates horizontally 
at a supersonic speed. The acoustic wavelength is L a = 2n/k a = 2ir/\/K 2 + I 2 = a w /T. Both 
vertical and horizontal wavelengths are larger. For readability, the wave and bottom amplitude 
is not to scale, and we have reduced the number of waves in the group from 104 to 10. Other 
than that, the angles are preserved. The configuration shown here corresponds to the conditions 
for maximum amplification of mode 1 (see below), with a vertical wavelength to water depth 
ratio of 0.75. For a water depth h = 4400 m depth, this corresponds to L = 7.7 km and a seismic 
frequency f s — 0.29 Hz. 



by figure [5J The coupling with the bottom motion selects a few resonant modes that 
dominate the solut ion, with an energy that grows linearly with the propagation distance 
( Hasselmannl [l~963 ) . For wavenumbers that allow a vertical propagation in the crust, we 



obtain body waves. These waves can be compression or shear (P or S) waves, and these 
two types occur in overlapping ranges of K. 

3.2. Acoustic noise in an unbounded ocean 

In order to illustrate the different types of solutions, it is interesting to evaluate the 
solution for an unbounded ocean, in which sound waves are radiated from the surface 
only. The velocity field and associated pressure fluctuations are 

fa = J_ f ™M K > / ) el [_„ +e(k , w , s y )] dKd fg (3 5) 

p w J -w; + igl 

P2 = f (K /,) el [- iz+e( k,k', s y)] dKd fg (3 6) 

J 1 - igl/uj 

where P2 has been obtained using the linearized version of eq. (|2.4[) . The measured 
pressure signal is the sum of the linear pressure p±, the second-order wave pressure P2 
given by eq. (|3.6p . and the Bernoulli correction P2,b given by 



P2,b( z )=P™ E ^(k, s ,k', s ',z)^ k ^; k ,e ie ^ k '' s ^'). (3.7) 



k,s,k' ,s' 



We note that p2,bot defined in eq. (|2.3ip is equal to P2,b{z — — h). 

We shall neglect g\l\/uj 2 s , which is bounded by the ratio between the deep water gravity 
and sound speeds, which is less than 0.1 for wave periods less than 180 s. We express the 
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velocity potential as a sum of propagating (acoustic, I real) and evanescent (acoustic- 
gravity, I imaginary) modes, 

02 = (j>2,p + 02, e- (3.8) 

We get the frequency spectrum of the propagating modes by integrating over the inner 
regions of the wavenumber space (labelled P+S, S and R in figure [1]), 

F p 2,ptfs)= I F p2 , S urf(K,/ 5 .)dK. (3.9) 

For this range of wavenumbers |fc — k'\ < K < uj s /a w , and using the relations lo s ~ Aitf 
and, (for small \ f — f'\), \k—k'\ ~ 2ir\f—f'\/C g ~ 87r 2 /|/ — /'|, we obtain an upper bound 
for the frequency difference |/ — /'| < g/{2-Ka w ) which is close to= 0.001 Hz. Typical 
ocean wave spectra have a relative frequency half- width tr /// that is between 0.03 for 
swells and and 0.07 for wind-seas (jrlasselmann et aZlll973l ). so that E(f) ~ E(f') is a 
good approximation for the interactions that drive long wavelength pressure fluctuations. 

The wave spectrum is thus broad enough for us to evaluate F P 2 <SUT { at K = using eq. 
(12.291) . and take it out of the integral in eq. (13.91) . The acoustic spectrum simplifies to 

W/ s ) = ^rPig 2 fsE 2 (f)i(f). (3.io) 

This is identical to the expression given bv lLlovdl (|l98lh . 



3.3. Gravity noise in an unbounded ocean 

The pressure associated with acoustic-gravity modes is the other part of the integral 
in (|3.9I) . for K > u) s /a w . The imaginary wave number I gives a vertical attenuation 
of the power spectrum by a factor e _2 ''K With that attenuation we may, for large 
enough depths, assume that only modes with K <C k contribute to the result, so that 
we may take F P 2, SU rf (K, / s ) ~ F p2 ,surf (K — 0, / s ), and take it out of the integrand. This 
approximation is valid only up to a maximum wave number AT max that is a small fraction 
of k, A'max = ek. For numerical applications we used e = 0.2. 
With this approximation we have, 

Fp2,e(f.,z) = ^ksurf ( K = °> fs)^ / K^AK 

K„ 



= F P 2, S urf (K = 0, / s )2tt f '" aX \l\e 2Wz d\l 
Jo 



^p 2 w g 2 f s [1 - e 2 **—] A 2 (/)/(/) (3.11) 



7T 

= 2? 

A previous investigation by Cox fc Jacobs! ( 1989h included an extra factor (1 + zK max ) 



in front of the exponential term e 2z - ft:max , because they neglected compressibility effects. 
That term , however, is negligible in the upper part of the water column, and their 
observations collected within 100 to 290 m of the surface in 4000 m depth, are thus not 
affected by this small compressibility correction. 

As shown in figure |B1 the oceanic pressure signals can be dominated by linear grav- 
ity waves down to depths of a few hundred meters. When looking at the double fre- 
quency band, linear waves may only dominate in the top 100 m. At these frequencies, 
the acoustic-gravity modes have the most important contribution between about 100 to 
500 m, provided that E 2 (f)I(f) is large enough. Propagating modes should dominate 
only beyond about 1000 m in the case of an unbounded ocean, or only 300 m, when 
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Figure 6. Example of (a) directional wave spectrum and (b) resulting profiles of the different 
contributions to the pressure fluctuations in the ocean, assuming infinite water depth. The ratio 
of double frequency to linear wave contributions depends on the amplitude of the waves and on 
the directional spectral shape, because all double frequency contributions are proportional to 
E 2 (f)I(f). This directional spectrum was estimated with a numerical wave model, and corre- 
sponds to the loudest noise event recorded at the ocean bottom seismometer H20, on 31 May 
2002 at 25 N, 136 W. This unusual spectrum has large wave energies in opposite directions, 
radiated from a Nor th Pacific storm and Hurricane Alma (This event is analyzed in detail by 
lObrebski et al.ll2012f ). 



accounting for the reverberation in a finite depth occan,assuming a typical tenfold am- 
plification for a sea floor with realistic elastic propertied. These depths will be reduced 
in the case of surface gravity waves with periods shorter than the 15 s swells example 
shown in figure O 

3.4. Atmospheric noise source: microbaroms 

The source of noise in the atmosphere can al so be derived with the sa me formalism, as 
an alternative to the Green functions used by Waxier fc Gilbert (|2006t . Indeed, we may 



consider the atmospheric motion to be irrotational, so that the equations of motion are 
identical in the atmosphere and in an unbounded ocean, with the only difference that 
the atmospheric density is p a and the atmospheric sound speed is a a - The second-order 
velocity potential takes the form, 



62,0 <x exp[i(K x x + K y y + l a z — w s i)] for z > 0, (3.12) 



with 




l a = J^--K\ (3.13) 



Because p w /p a — 1000, the air motion has only a small 0(p w /p a ) local influence on 
the water motion, so that the solutions derived earlier for the water motion remain valid 
in the presence of air. The air motion, with a velocity potential <p a also obeying eq. (|2.13|) 

f This amplification depends not only on the impedance ratio of the water and crust, which 
defines the amplification coefficients Cj derived below, but also on the sei smic attenuation coef - 
ficient Q, which is discussed in section 4. Realistic calculations following Ard huin et al\ (|201lh 
typically give a factor 10 to 20 amplification of the sound in the water column due to the bottom 
elasticity. 
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is fully determined from the water motion via the kinematic boundary conditions on the 
air and water-sides of the interface ()2.7|) . 



d(j) a dcf) d 2 (4> a - <f)) 

— — — ~ v (fa - 4>) • vc - C — ^ at z = °- ( 3 - 14 ) 

oz oz o i z 
From the first order potential in the air (e.g. lWaxler fe Gilbert 2006 ) 

<pi=J2 is-^i, k e~ fc V (k-x ~ s,T) (3.15) 
we obtain the second order potential, 

d -^- d 4l + £ ^a(k, S ,k^ S ^z)^ k Z 1 Ve ie ^'-^ (3.16) 
and a new coupling coefficient 



D za (k, s, k', s', z) = -^- [kk 1 + k • k') 

(7 



(3.17) 



We note that for k' = k, D za = 0, so that the long-wavelength motion with K <^k 
simplifies to 

d(j>2a 9(j>2 



dz dz 



at z = 0. (3.18) 



consistent with the result given by iPosmentierl (|1967l ) for the interaction of monochro- 
matic wave trains, and in disagreement with a factor 8 correction proposed bv lArendt fe Frittd 
(|2000h . 

This gives a pressure spectrum for the propagating atmospheric waves, 

F p2 , ap (fs) = I 4fi^ P 2, S urf(K, /.)dK = R(K m ^)^plg 2 f s E 2 (f)I(f). 

(3.19) 

with the non-dimensional factor 

2 rK max 1,21 2 rK max 1,21 

R(K max ) = 2^ / ^KdK == 2=f / y XdX. (3.20) 

In order to avoid the singularity for l a = 0, and atmospheric ducting effects not repre- 
sented here, we take K max = oj s / (2a a ) which restricts the acoustic propagation directions 
to be within 30 degrees from the vertical. In that case we have R(K m ? x) ~ 0.54 instead 
of fitf inax) = 0.25 with the vertical propagation approximation of IWaxler fc Gilbert! 
(2006), which replaces the l 2 /I 2 , factor in the integral by its value a 2 /a 2 u for K = 0. 
Other than th at, our expression is con sistent with their low Mach number asymptote, 
i.e. a/k <C a a (|Waxler fc Gilbert! |2006[ eq. 61). The present theory also allows the esti- 
mation of the evanescent wave components given by wavenumbers K > u> s /a a . 



4. Noise in a finite depth ocean 

For large depths compared to the OSGW wavelength, kh ^> 1, the finite depth has 
little effect on the evanescent modes except for a doubling of the motion amplitude 
near the bottom, as the vertical profiles of the form exp(Kz) are replaced by cosh(Kz + 
K K)j cosh(Kh) . This is similar to the finite depth effect on linear wave motions. However, 
the propagating modes radiated by the surface will now undergo multiple reflections at 
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the bottom and sea surface, as shown in figure [5] The oceanic acoustic field is tightly 
coupled to elastic waves in the crust through these reflections. 

One of the greatest complications induced by the presence of a bottom is the hetero- 
geneity of the sediment and rock layers below the water column. The natural layering 
of the crust has a strong influence on the sound reflection and the nature of the seis- 
mic modes (e.g. Latham fc Suttonl 119661 Abramovicil 19681 ). T hese effects will no t be 



considered here, and we follow exactly the theoretical setting of Hasselmann ( 19631 ) . 



4.1. Elastic wave theory 

For simplicity, we assume here that the ocean bottom is a uniform and semi-infinite solid, 
with constant density p s and compression and shear wave velocities a c and j3 given from 
the Lame coefficients A and /i, 



■2 A + 2/i 2 _ p 



<4 = -^Z, F = ^. (4.1a,6) 

Ps Ps 

Assuming that the crust is a Poisson solid, fj, = A and a c = y/3/3. 

The free wave problem in these con ditions was solved by IStonelevI (|l926h . Here we 



consider the forced problem treated by lHasselmann ( 19631) with a forcing by a pressure 



field P2,surf at the sea surface, but now generalized to an additional wave-induced bottom 
pressure p 2 ,bot- 

The equations of motion in the water column are unchanged from the previous section, 
but they are now coupled to the elastic motions of the crust. Crustal motions can be 
separated in an irrotational part with a velocity potential 4>c and a rotational part with 
a stream function ip, both solutions to Laplace's equation. With a wave source at the 
surface, a horizontal propagating wave of phase O = Kx — uj s t, implies that cf> and ip are 
either decaying or propagating downwards. They must therefore take the following form, 

<p c = Ae rn( - z+h ^e ie , ip = Be n( - z+ Ve i@ . (4.2a,&) 

Both vertical wavenumbcrs m and n are given by the Fourier transform of eq. (|2.13p . 



a 




m= Jk 2 -^f, n = JK 2 -^ (4.3o,&) 



where the sound speed in water has been replaced by the compression and shear velocities. 
For K > ui s /P s , m and n are real and both compression and shear waves are evanescent. 
For uj s /a c < K < (3 S the compression wave is evanescent but there is a shear (S) wave 
that propagates through the crust. 

The constants A and B have dimensions of m 2 /s and are determined by the boundary 
conditions at the ocean bottom. 

Horizontal and vertical ground displacements are given by the real parts of 

£r = (KAc m ( z+h *> + inBe n{z+h ^ e ie /oj s (4.4) 
& = (-imAe m{z+h) + KBe n{z+h) ^ c ie /uj s . (4.5) 
Hooke's law of elasticity gives 

T zz — A I — h — I + 2^—, t xz = p — + — 4.6a,6 
\ox oz J az \ az ox J 

The zero tangential stress on the ocean bottom t xz (z = — h) = yields the following 
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2\Km 



B 



17 



(4.7) 



n 2 + K 2 ' 

Thus, in addition to the unknown water-side amplitudes C and D of the velocity 
potential at the bottom, we have one more unknown, the compression wave amplitude 
A on the solid side. 

The three equations that relate C, D and A are: the combined kinematic and dynamic 
boundary condition (|2.15|) . and the bottom continuity of normal velocity 



dz 



d(j) c dtfj 

— h — at z= -h 

oz ox 



and normal stress, 



T zz {~h) = p(-h) = p n 



dt 



P2,bot- 



(4.8) 



(4.9) 



For waves in intermed iate or shallow wat er, i.e. kh < n, the Bernoulli term p2,bot that 
was not considered by lHasselmann ()l963f ) should be included. We thus obtain the linear 
system of equations 



■gii) 



e Uh C 



qA - UC 
rA - iuj s p w C 



gil) e" 

+ilD 
iu s p w D 



«^D = ^p 2 , surf (K,/ s ) 

Pw 

= 



-p2,bot(K,/s) 



with 



i 

UJ S 

i 



2K 2 p 2 

4:/3 4 K 2 mn 
~ lo 2 s -2K 2 I3 2 

p s m 2 a 2 c + \k 2 



2o2\ 



+ {luI - 2K 2 (3 
K 2 mn 

+ 4^i- 



K 2 



(4.10) 

(4.11) 
(4.12) 



(4.13) 



(4.14) 



Since we are in the range where I < uj s /a w we may neglect again g\l\/uj 2 , which is less 
than 0.1 for OSGW periods less than 180 s. We rewrite these equations in matrix form, 

M[A, C, D] T = [-p 2 ,8urf(K,/ a ), 0,-p 2 , bot (K,/ s )] T 



with 



M 



ip w uj s e llh ip w uj s e uh 



(4.15) 
(4.16) 



q —il il 

r -\L0 s p w ~iuJ s p v 

The general solution of eq. (|4.15[) is the sum of one particular forced solution and the 
general solution of the homogeneous system, without the right hand side forcing, i.e. the 
free waves. This combination of free and forced waves is fully determined by the initial 
conditions. 

The forced solution is readily expressed using the determinant of the system 

Uh. 



det(M) = -p. u 



(le ilh r + ui s p w e ilh q~le [lh r - w s p w e"- a ) 



2p v 



2K 2 fi 2 ) 



[iPs 



cos(//i) U[3 A K 2 mn - (lu 2 - 2K 2 [3 2 ) 2 



(4.17) 

— p w m sin(//i)o>g| , 
(4.18) 
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in the form 

A = 21 Pw oj s det(M) ( 4 ' 19fl ) 

= -i (lr + ga; s p tu )p2, S urf (K, / s ) + ge~'^p2 : bot.(K, / s ) 

det(M) 1 ' ' 

D = i (qu sPw - Ir) ^ = Zr > P^|lMz gg^Wj^ /') . (4 .l 9c ) 

As detailed below, this determinant vanishes for the pairs (u> s ,K) that fall on the dis- 
persion relation of Rayleigh waves. But the solution fo r random waves ca n always be 
obtained by integrating across this singularity, following Hasselmannl (1962). 
For kh^> 1 we may neglect p2,bot and, 

c _ P2,surf(K,/ s ) ir' - q' (4 20a) 

p w tu s 2r' cos(lh) — 2q' sin(Zft) 

D= P2,snrf(KJ s ) ir'+q' 

p w u! s 2r' cos(lh) — 2q' sm(lh) 



with 



p w Lo s q, r = ilr. (4.21 a, b) 



4.2. Acoustic noise in a bounded ocean 

Taking typical values of the water and crust density and sound speeds gives r' /q' > 
70 for the free modes that are significantly generated by the waves (i.e. Cj > 0.1, as 
defined below by eq. I4.34|) . We can thus consider that q' /r' <C 1, which gives C ~ D ~ 
P2,surf(K, f s )/[2p w oj s (cos(lh) + q'/r')]. The velocity potential and pressure in the water 
are given by, 

, . R fv , > cos[l{z + h) -q'/r'] ie cos[Z(z + fc)] - q'/r' ie 

02 ^ ip 2 ,surf K, / s ) — — — e , p 2 ^ ip 2 ,surf K, f s )— — —— — e . 

p w u} s {cos(lh) — q'/r) (cos(m) — q'/r') 

(4.22o,6) 

The small but finite factor q'/r' ensures that the solution remains finite as a small 
fraction of the acoustic energy is radiated into the crust, otherwise the acoustic energy 
would accumulate in the water column. The pressure oscillations are thus maximum 
for resonant frequencies such that the ratio of the water depth and vertical wavelength 
lh/{2%) is 1/4, 3/4, 5/4 ... We note that the vertical wavelength 2ir/l is always greater 
than the acoustic wavelength 2n / \JK 2 + I 2 — f s a. As a result, the resonant frequencies 
are shifted to higher values compared to the vertical resonant condition that is given by 
f s ha w = 1/4 ... 

We may now integrate the pressure spectrum for all acoustic wavenumbers to find, 
again, the frequency spectrum, 



F p2iap (f s )=27Tp 2 w g 2 f s E\f)I(f) 



r cc . s (; z + ^)^ g y r ' -' 2 

cos(lh) ~ q'/r' 



KdK. (4.23) 



We now have a depth dependence of th e sound spec trum. At the surface it is equal to 
the widely used unbounded ocean value ()Llovd 1981 ) , but in the water column it can be 
strongly amplified at depths where cos(//i) approaches 0, which includes the near-bottom 
region where all the resonant modes have an anti-node. 
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Figure 7. Wavenumbers of seismic and acoustic modes for a fixed frequency, which are the 
solutions of eq. (|4.24|) . The right hand side of that equation is abbreviated as RHS. Here we 
used h = 5000 m, a w = 1500 m/s, /3 = 3200m/s, p w = 1000 kg/m 3 and p s = 2500 kg/m 3 . This 
graph corresponds to a seismic frequency f a = 0.263 Hz, for which f„h/a w — 0.88. Here the 
determinant has three real roots, one acoustic mode K a for which I = (horizontal propagation) 
and two seismic modes Ko and K\. The number of these roots increases with the frequency. 
One new solution appears every time 27r/ s becomes larger than cj°j defined by eq. (|4.29[) . For 
reference we also give the wavenumber for which the right hand side of eq. (|4.24|l is zero, which 
is the Rayleigh wavenumber Kr in the absence of the ocean layer (h — 0). 



4.3. Rayleigh waves 

In order to simplify the algebra, we consider in this section waves in deep water and 
neglect bott o m for cing. As a result, this section brings no new results compared to 
Hasselmannl ( 19631 ). but the properties are discussed in more detail and more explicit 



expressions are presented that will be used later in the derivation of new solutions for 
other types of seismic waves. 

For a fixed frequency uj s , there is at least one wavenumber K for which det(M) = 0. 
This condition defines the dispersion relation of the Rayleigh modes ( Stonelev 1926), 



, . lp s 4p 4 K 2 mn - (u 2 - 2K 2 /3 2 ) 

tan (Ih) = x — V ( 4 - 24 ) 

mp w u* 

with the fundamental mode corresponding to the largest K value. Figure [7] illustrates 
the family of modes for a give frequency. 

The phase speeds of the Rayleigh modes vary continuously from the shear wave velocity 
f3, in the limit n = where the shear waves transition from propagating to evanescent, 
to the sound speed in water a w , in the limit I — where the acoustic modes become 
evanescent in the water (figure [8] a). This variation of the phase speed has an inflexion 
point close to the phase speed of Rayleigh waves without the water layer, corresponding 
to a maximum in group speed (figure E]b). 

We may now use eq. (|4.19p to obtain the ground displacement amplitude 



5 = £z{-h) = G(K Xl K y ,uj s )p2. S u 1 -f{K x ,K y ,uj s ) 



(4.25) 
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1000 



Figure 8. Phase (top) and group (bottom) speeds of Rayleigh waves as a function of the 
dimensionless water depth. Calculations used a w = 1470 m s -1 , /3 = 3000m s~\ p w = 1000 and 
p s = 2500 kg m -3 . Values without the water layer are also indicated with the blue line. In that 
case Rayleigh waves are not dispersive. 



as a function of the amplitude of the sea surface pressure, with the transfer function 

2ip w lm(K 2 -n 2 ) U26) 
(n 2 + K 2 ) det(M) ' ^ ' ' 

An example of this transfer function is shown in figure GO in which the dispersion 
relation of the Rayleigh modes appear as the narrow red bands around the singularities 
of G. 

To evaluate the complete family of Rayleigh wave solutions with dispersion relations 
uj s j (K) , where j is the mode number, we need to examine the nature of these singularities. 
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Figure 9. Magnitude of the transfer function G(K,f a ) illustrating the singularities along the 
dispersion relations of the free Rayleigh modes with wavenumbers K — K n (f a ), that correspond 
to det(M) = 0. The oblique dashed lines correspond to phase speeds equal to ct c , /3, and a w , and 
separate the four domains body waves (P+S), mixed body and evanescent waves (S), Rayleigh 
waves (R) and acoustic-gravity modes (A-G). 



We re-write the determinant of the system as 
2ilp s p w cos(lh) 



det(M) = 



{ (4/3 2 /c 2 ) (mn + u% - p 2 K 2 ) - uj 4 s (1 + tan(lh)p w m/p s l) } 



(wf - 2K 2 /3 2 ) 

(4.27) 

For most frequencies, the singularities are simple, allowing a Taylor expansion of G in 
this form 

G'{K 3 (lu s )) 



G{K,u s 



However, for each mode j, there is a critical frequency 



~c to „<-A w „> 
m 

2«2 



■ J7T 



(4.28) 



(4.29) 



for which both t&n(lh)p w m/ p s l = 1, and w 2 — K j3 = 0. These critical frequencies 
exist for modes j > 0. These frequencies are those for which new Rayleigh modes 



9, J 



appear, in a way similar to the Love waves discussed by Aki & Richards (2002, figure 
7.3). The singularities at (K,uj s ) = (ui®j//3,ujgj) are not simple and for these we have 

G(uj s , k) ~ G'iKAu.y/yjK* - Kjfa). 

Solutions for the vertical displacement at the top of the crust have a spectral amplitude 
which is linearly related to the equivalent surface pressure amplitude, 

S(K x ,K y ,uj s ) = G(K X , Ky,u) s )p 2<aui {(K x , K v ,u a ). (4.30) 

For the simple singularities of G, we may write G(K,uj s ) — G'(K, uj s j)/(uj 2 — w^-) + 
0(uj s — uj s j). Taking initial conditions S(t = 0) = and 85/dt = gives the full solution 
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( Hasselmannlll962l eq. 3.2) of the j th Rayleigh mode response 



5{K x ,K y ,ui s )=G'{K,u> s A 



-iui s t 



2uj, 



(4.31) 

as a function of K, where (J s j(K) is frequency of the j th Rayleigh mode. For a forcing 
that varies slowly on the scale of the seismic period 2it/uj s , and provided that the forcing 
spectrum is continuous in spect ral space, th i s give s the rate of change of the ground 
displacement spectrum given by lHasselmann ( 1963f ). 



dF s (K x ,K y 
dt 



S DF (K x ,K y ) 



tt\G 



'i 2 



2w 2 ■ 



--Fp2,surf (K x , K y , U s ). 



(4.32) 



The other discrete singularities, at (K, oj s ) = -//3, w° ■), are associated with 'con- 
ical' or 'head' waves (|Aki fc Richards! 120021 ). for which the vertical wavenumber n = 0, 
and that propagate along the ocean-crust interface. These horizontally propagating shear 
waves are generated with an evanescent compression wave. The singularity is integrable 
over the two spectral dimensions lo s and K . 

Although it looks like only the resonant forcing contributes to the solution, it is in 



fact the near-resonant forcing (uj s 



J s,3 



) that builds up the seismic noise, because the 



exact resonant terms have a zero measure in spectral space. Indeed, a purely resonant 
forcing would give an amplitude that increases linearly with time, and an energy that 
increases like t 2 . The linear growth of energy in time can be interpreted as an effect of the 
narrowing with time of the frequency bandwidth in which th e interaction is sig nificant. 
This is a general property of wave- wave interactions (see also lHasselmannl ll966). 

We express the source of seismic noise with the rate of increase of the variance of 5 
per unit of propagation distance. It is 



Sdf(w s ) 



K(s)S DF (K x ,K y ) _ ^ 2 f s c 2 3 



U 2 



Fp2,surf(K x ,K y , UJ S ) 



(4.33) 



where U is the group speed of the seismic waves, and Cj is a dimensionless coefficient 
that depends on uj s h/a w and the seismic mode index j, shown in figure [TU1 



U 2 2ttuj s 



2uj 2 



(4.34) 



We note that a missing 2ir in eq. (5) of I Ardhuin et all (|2011l ) has been corrected here. 

A very rough simplification can be obtained by taking U and K independent of j and 
h. Then Sdf(/s) can be taken as a sum of all Rayleigh modes in the form 



5'df(/s 



(4.35) 



Values of Cj are obtained from eq. (|4.34l) , and shown in figure [TUJ 

We may propagate these sources of seismic waves in a vertically symmetric earth model, 
neglecting all three-dimensional propagation effects, and parameterizing seismic wave 
scattering and dissipation with a uniform quality factor Q. Under these assumptions, the 
spectral density of the vertical ground displacement at z = — h and at the longitude Xq 



double-frequency noise generation by surface gravity waves waves 



23 




Figure 10. Dimensionless coefficients Cj that amplify the wave-induced pressure into ground 
displacement. The maxima of Cj correspond to quarter-wavelength resonance typical of organ 
pipes, except that here the sound waves propagate obliquely in the water column, which is why 
the maxima are at values of f s h/a m which are not exactly at 1/4, 3/4..., but shifted to higher 
frequencies by a factor 1/ \f (K 2 + l 2 ). The amplitudes of the peaks depend on the impedance ratio 
of the sea water and crust. Hence the peak amplitude increases with p 3 /3/(p w a w ). For example, 
f3 — 2800 m s _1 , gives a maximum of 0.88 for Co instead of 1.03 here with /3 = 3000 m s _1 . 



and latitude <po 



r/2 



-7T/2 



F s (\o,<f>o,f s )= I' I |^4%e- 2 ^^/(^)(i?| S in0' o dA' o d^ o ) (4.36) 



Re sin A 



with Re the Earth radius, and U the seismic group velocity. The term {R 2 E sin (^dA^d^) 
is the Earth surface area element. The denominator (Re sin A) is the geometrical spread- 
ing fa ctor for wave energy that follows geodesies on the sphere (e.g. iKanamori fc Given! 
1981 ). replacing the distance (i?eA) used in flat Earth models (e.g. Hasselmann 19631) . 



The Rayleigh waves thus generated propagate away like free modes, those that exist 
without the local forcing. For these free modes with a monochromatic ground displace- 
ment S(x, y, t) the surface pressure is constant so that C = —D and the velocity potential 
and pressure in the water take the form, 

fa(z)= —. sin(lz)5(x,y,t) (4.37a,6) 
/ cos(lh) 

uj 2 

P2{z) = Ww-, tttt sm(lz)8(x, y, t). (4.37a,6) 

/ cos(ln) 

For varying water depths, and on land, one may assume that t he seismic energy i s 
propagated along a ray and apply the refraction coefficient given by lHasselmann (1963). 



4.4. Seismic P and S waves 

Unlike the Rayleigh waves that grow resonantly in time owing to the trapping in the 
ocean/crust waveguide, the P and S waves, for which the vertical wavenumbers m or 
n are complex, radiate into the earth's interior, and their energy level is given directly 
by its value at the source. In particular, for K < u! s /a c we have propagating P waves 
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Figure 11. Non-dimensionless coefficients cp and cs that amplify the wave-induced pressure 
into ground displacement associated with P and 5* waves. 



with a velocity potential amplitude A given by eq. (I4.19[) . In the particular case K = 0, 
which correspond to standing waves, we have only have P waves, no S waves or Rayleigh 
waves, and these propagate exactly along the vertical axis. The frequency spectrum of 
vertical ground displacements at (z = —h) can be evaluated directly with (|4.30[) because 
G has no singularity in this range of wavenumbers, 



Fs, P (fs) = f s E 2 (f)I(f) 



2 2 



with a non-dimensional coefficient cp, 



•A 



4l 2 m 2 P 2 s (3 4 



-KAK. 



(4.38) 



(4.39) 



w s 2 det 2 (Af)' 

A similar expression can be written for 5* waves, and both are illustrated in figure [TT1 



Fs.sifs 



f s E 2 (f)I(f)^c 2 s . 

PsPs 



(4.40) 



However, in the range of wavenumbers where S waves exist, k < lo s / 1 (3, there can also 
be evanescent P waves, and the system can approach the singularity for lo s — w SJ - and 
k = oj s //3. We evaluated numerically the coefficient 



4 = 2n 



4Pm 2 k 2 p 2 p 4 



uj 2 s (n 2 +fc 2 )det/(M) 



-KAK. 



(4.41) 



Due to the typically th ree times stronger attenuation of S waves com pared to P waves in 
the Earth mantle (e.g. Anderson fc Hartlll978t [Pasvanos et al. 20091 ). we will now focus 
on P waves only, which should dominate in the far field of the noise source. 

For the estimation of the spectrum recorded outside of a source area, it is more con- 
venient to express the local seismic source as a function of the horizontal propagation 
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Figure 12. Dimensionless coefficients cp ltp that amplify the wave-induced pressure into ground 
displacement. The maximum for a zero take-off angle corresponds to vertically propagating 
compression waves, and the compression waves that propagate along the crust have a vanishing 
amplitude. 



angle 8, and the vertical take-off angle ip. For P waves, this gives, 
F s ,p(f a) 9,<p) = f s E 2 (f)I(f)^-c P ^sm^ 

PsPs 

with the non-dimensional coefficient cp^ defined by 

a _ 4i 2 m 2 ^/? s 4 dK 



-p,<p 



uj s a c det 2 (M) dip ' 



(4.42) 



(4.43) 



which is the normalized source per unit solid angle J7, so that the average over the half 
space of downward directions fi - is 



2 f 2 * r' 2 2 f 2 

c P = I / c P ^ sin ipdipdB = / c P tp dQ, 

Jo Jo Jn- 



(4.44) 



as defined by eq. (I4.39|) . 

It is noteworthy that the distribution of the P-wave energy with the take-off angle is 
very clos e to the one given by a small disk pu shing at the t op of a uniform half space, as 
given bv lMiller fc Pursev (1955) and used bv lVinnikl ( 1973 ). although it also vari es with 
the no n-dimensional water depth f s h/a w . The only missing item in the work by I Vinnik 
( 19731) is the very strong amplification of the motion for resonant frequencies associated 
with the water layer. Due to the large impedance contrast at the water-crust interface the 
relative amplification of P waves is one order of magnitude stronger than for Rayleigh 
waves. We thus expect a much tighter correspondence of the strong seismic noise sources 
with the water depths that correspond to a maximum amplification. 

4.5. Observable P wave spectra 

We will now finish our analysis of these body waves by estimating the ground motion due 
to P and Rayleigh waves as a function of the distance from the source, an application 
of practical interest. For a seismic station or hydrophone in the ocean, the incoming 
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energy at the receiver will amplified by multiple bottom and surface reflections in the 
water column. However, this P-wave signal from remote sources is likely to be dwarfed 
by locally generated noise. 

For a land-based station, we may assume that P waves arrive directly from the source 
area. There are many other seismic wave phases that have undergone multiple reflections 
at the surface, these are called PP, PPP or at any inner boundary of the Earth , 



like PkP phases that have reflected off the mantle-core interface (|Aki fc Richardsll2002l ). 
These can be treated exactly like the direct P phase. 

On arrival at the receiver, these waves are totally reflected at the Earth surface at 
z = 0, which doubles the ground motion, so that the ground displacement is given by the 
integral of four times the incoming spectral densities over the directions 9 and (p, which 
can be replaced by an integral over the source positions (As, (j>s)- The transformation from 
the ray parameters (#, tp) to the geographical coordinate can be obta ined approximately 
for any type of seismic wave using travel time tables (e.g. lSnokdl2009l ). which also provide 
the travel time r. 

We may now express the ground displacement due to P waves at the observing station 
of coordinates (Ao, <Po)i as a function of the same quantity at the location of sources, as 
given by eq. (|4.42l) 

*i,p(Wo,/.)= / ^ P (\^f s ,9^)e- 2 *f^ s ^IQ<m, (4.45) 

with dtt an element of solid angle that corresponds to the ensemble of rays arriving from 
an Earth surface element around the sources located at (As, <f>s)- This sum may also be 
transformed as an integral over the ocean surface by properly mapping £1 to (\s,<j>s)- 
The elementary solid angle df2 is zero for the so-called shadow zones, the regions for 
which there is no P-wave ray that connects to the observing station. For a single phase 
of seismic waves, this ensemble of rays has a half-banana shape, as illustrated on figure 

From our calculations, we expect that P waves will dominate the signal at large dis- 
tances from the source. The exact location where P-wave levels overtakes Rayleigh-wave 
levels depends on the relative attenuation of the two types of waves. With a realistic 
Q = 2000 for the P waves, and Q = 400 for the Rayleigh waves, figure [T4l shows that it 
occurs at an epicentric a ngle of 40°, wh ich is a distance of 4400 km, consistent with the 
observations reported bv lVinnikl (|l973l) using Kazakhstan array data. 



5. Conclusions 

We have shown how the same physical process, the interaction of ocean surface gravity 
wave (OSGW) trains, can produce a wide variety of noises, in the atmosphere, ocean 
and Earth's crust, that can be classified according to their horizontal phase speed. The 
slowest noises in the ocean are acoustic-gravity waves that dominate pressure records 
at depths less than about one tenth of the acoustic wavelength. These acoustic-gravity 
waves cannot exist in the absence of OSGWs and are thus confined to the region of active 
wave forcing. Intermediate phase speeds correspond to Rayleigh waves that contain most 
of the energy of the seismic modes for distances less than about 4000 km from the source. 

We corrected previously published asymptotic behaviour for very long period noise 
(T > 30 s). In particular we find that the sources of this long period noise are attenuated 
on the continental shelves, consistent with previous studies of forced gravity wave motion. 
This finding supports a spatial distribution of these sources outside of the continental 
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Figure 13. Schematic of rays for seismic P waves radiated from a point source (red dashed 
lines) with a directional distribution cp tip (solid red line), and ensemble of rays received at a 
given station from an extended noise source (yellow half-banana). 
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Figure 14. Estimates of the rms vertical ground displacement associated with Rayleigh 
or P waves, as a function of the epicentric angle A, for a source of intensity 
/ Jp2,sm-f(K = 0, /)d/ = 4.2 x 10 4 hPa 2 m 2 over a 330 by 330 km square, assuming an at- 
tenuation factor Q = 2000 for th e P waves ( Pasyanos et al. 20(3), with travel times given by 
the akl35 reference Earth model (Snoke 2009). 
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shelve s, on the shelf breaks, which is co nsistent with data from Rhie fc Romanowicd 
(2006) or in deeper water, as reported by iNishida fe Fukaol (|2006l ). 

The common source of all these noises should allow a verification of the source magni- 
tude for seismic waves by near-surface measurements of pressure which is dominated by 
acoustic-gravity modes. In particular, the direct modelling of the acoustic-gravity modes 
can be compared to pressure measurements i n depths less than a few hundreds of me- 



ters. Unlike the analysis of seismic noise (e.g. I Ardhuin et aLll2011l) . which suffers from 



poorly known seismic propagation and attenuation factors, th e acoustic-gravity wave 
attenuation over the water column can be predicted accurately (jHerbers fe Guzalll994T) 
and thus pressure measurements in the upper ocean may provide a more quantitative 
verification o f numerically modeled direc t ional surface wave properties. In particular, as 
proposed by ICooper fc Longuet-Higgins (1951), pressure measurements may provide a 
precise estimate of coastal reflection or wave scattering by currents, sea ice or other ef- 
fects. Noise records from land-based or seafloor-mounted seismometers are more ambigu- 
ous because they integrat e sou rces over a large area. Al so, as discussed by IHasselmann 
( 19631 ) , Abramovici ( 1968 ) , and Latham fc Sutton ( 1966 ) , the variations in water depths 
and horizontal and vertical variations of properties in the Earth's crust can significantly 
modify noise properties in both the water column and the crust. 

F.A. is funded by ERC grant #240009 "IOWAGA" with additional support from the 
U.S. National Ocean Partnership Program, under grant N00014-10-1-0383. T. H. C. H. is 
supported by the U.S. Office of Naval Research Littoral Geosciences & Optics Program 
and the U.S. National Oceanographic Partnership Program (NOPP). Discussions with 
S. Webb and L. Marie on the theoretical aspects are gratefully acknowledged. 



Appendix A. Coupling coefficients 

Using the coupling coefficient D given by IHasselmann ( 1962 , eq. 4.3) for the velocity 
potentials, our coupling coefficient for the elevation amplitudes is 



D z (k, s,k',s' 



g 2 J(k, S ,k', S ') 
isas'a 1 (sa + s'a') 



k-k' 



0.5 



Vfc 2 \ 



(sa + s'a') Vcosh 2 (fc'/i) ' cosh 2 (fc/i)/ 

(Al) 

In the bottom pressure, the additional term arising from the orbital velocity has a 
coupling coefficient 

D p b (k, s, k', s', z) = 

2 kk' sinh[fc(z + h)] sinh[fc'(z + h)] - k • k' cosh[£;(z + h)] cosh[fc'(z + h)] 



2sas'a' cosh(fc/i) cosh(k'h) 



(A 2) 



The relationship with the coupling coefficient C given by iHerbers fc Guza (Il99ll . their 
eq. 4) for the bottom pressure, expressed in meters of water, is given by solving eq. (|2.15l) 
for </> 2 , and then rewriting Bernoulli's equation (I2.4[) . as 



EL 

Pw 



dt 



|V0x| 



dz 



(A3) 
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This gives, for z = —h, 

c _ D^sa + s'a'f D pb (z = -h) 

g[gKtarLh{Kh)-(sa + s'cr>) 2 } g [ ' 
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Symbol meaning where 



1 


and 2 


indices denoting first and second order motions 




a w 


sound speed in water 




a c 


compression wave speed in crust 




a a 


sound speed in atmosphere 




P 


shear wave speed in crust 




A 


angular distance 




5 


vertical ground displacement at the top of the crust, 5 = £ z (z = —h) 




c 


elevation of the sea surface 


A 


and fi 


Lame elastic coefficients for the crust 


Ao 


and <j)o 


longitude and latitude of the observation location 


As 


and (f>s 


longitude and latitude of a noise source 




P 


perturbation of density 




Pw 


mean water density 




a 


radian frequency of surface gravity waves a = 2nf 




T 


stress tensor 




e 


phase function of the seismic or acoustic waves 




<t>a, 4>c 


velocity potentials in the water, atmosphere, crust 




1> 


stream function in the crust 




<P 


take-off angle for seismic body waves 






displacement of particles 






radian frequency of noise u> s = 2nf 3 




<3 


critical values of oj s for which new modes appear when uj s increases 



Table 2. Table of arab and greek symbols 



Appendix B. Definition of symbols 
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eq 


g2 




eq 


(E3 




eq. 


12.11) 


eq. 





Symbol meaning where 

A, B amplitudes of <j>c, ip 

C, D amplitudes of up- and downward propagating components of 

C g Group speed of surface gravity waves 

D z coupling coefficient for the surface elevations 

a and a amplitude of surface gravity waves 

Cj non-dimensional amplification factor for Rayleigh mode number j 

cp and cs non-dimensional amplification factor for P or S waves 

/ and /' frequency of surface gravity waves 

f„ acoustic or seismic frequency 

_F P 2,surf spectral density of p 2 ,surf 

F5 spectral density of 8 

G surface pressure to bottom vertical displacement transfer function 

g apparent gravity acceleration 

h water depth 

/ directional integral of the wave spectrum eq. (|2.27[1 

j Rayleigh wave mode number, counted from 

k and k' vector wavenumbers of surface gravity waves 

K horizontal vector wavenumbers of acoustic or seismic waves 

I, l a , m and n vertical wavenumbers for <j>2,4>2,a , 4>c and ip 

p pressure 

P2,surf wave-induced forcing at the sea surface eq. (|2.16|) 

j?2,B Bernoulli head pressure eq. (|3.7[) 

P2,bot wave-induced forcing at the bottom, pb.bot =P2,b(z = —h) eq. (|2.31|) 

M matrix of the linear system of equations 

Q seismic quality factor (i.e. damping coefficient) 

q and r coefficients eq. ()4.14|) 

ijfffnjx) non-dimensional coefficient in microbarom source eq. 1)3.201) 

Re radius of the Earth 

Sdf Seismic source of Rayleigh waves 

s and s' sign variables equal to -f or 1 

u and w horizontal and vertical velocity components 

V group speed of seismic waves 
Zi k complex amplitude of the linear surface elevation component (k, s) 



Table 3. Table of roman notations 



32 



F. Ardhuin and T. H. C. Herbers 



REFERENCES 

Abramovici, Flavian 1968 Diagnostic diagrams and transfer functions for oceanic wave-guides. 

Bull. Seism. Soc. Am. 58 (1), 426-456. 
Aki, K. & Richards, P. G. 2002 Quantitative seismology. Sausalito, California: University 

Science Books, second edition. 
Anderson, D. L. & Hart, R. S. 1978 Q of the Earth. J. Geophys. Res. 83 (B12), 5869-5882. 
Ardhuin, F., Chapron, B. & Collard, F. 2009 Observation of swell dissipation across oceans. 

Geophys. Res. Lett. 36, L06607. 
Ardhuin, F., Stutzmann, E., Schimmel, M. & Mangeney, A. 2011 Ocean wave sources of 

seismic noise. J. Geophys. Res. 116, C09004. 
Arendt, S. & Fritts, D. C. 2000 Acoustic radiation by ocean surface waves. J. Fluid Mech. 

415, 1-21. 

Arrowsmith, S. J., Johnson, J. B., Drob, D. P., & Hedlin, M. A. H. 2010 The seismoa- 
coustic wavefield: a new paradigm in studying geophysical phenomena. Rev. of Geophys. 
48, RG4003. 

Bonnefoy-Claudet, S., Cotton, F. & Bard, P.-Y. 2006 The nature of noise wavefield and 

its applications for site effects studies a literature review. Earth Sci. Rev. 79, 205-227. 
Cooper, R. I. B. & Longuet-Higgins, Michael S. 1951 An experimental study of the 

pressure variations in standing water waves. Proc. Roy. Soc. Lond. A 206, 426-435. 
Cox, C. S. & Jacobs, D. C. 1989 Cartesian diver observations of double frequency pressure 

fluctuations in the upper levels of the ocean. Geophys. Res. Lett. 16 (8), 807-810. 
Duennebier, F. K., Lukas, R., Nosal, E.-Ma., Aucan, J. & Weller, R. A. 2012 Wind, 

waves, and acoustic background levels at station aloha. J. Geophys. Res. 117, C03017. 
Farrell, W. E. & Munk, W. 2008 What do deep sea pressure fluctuations tell about short 

surface waves? Geophys. Res. Lett. 35 (7), L19605. 
Farrell, W. E. & Munk, W. 2010 Booms and busts in the deep. J. Phys. Oceanogr. 40 (9), 

2159-2169. 

Fukao, Y., Nishida, K. & Kobayashi, N. 2010 Seafloor topography, ocean infragravity waves, 
and background Love and Rayleigh waves. J. Geophys. Res. 112 (10), B04302. 

Hasselmann, K. 1962 On the non-linear energy transfer in a gravity wave spectrum, part 1: 
general theory. J. Fluid Mech. 12, 481-501. 

Hasselmann, K. 1963 A statistical analysis of the generation of microseisms. Rev. of Geophys. 
1 (2), 177-210. 

Hasselmann, K. 1966 Feynman diagrams and interaction rules of wave- wave scattering pro- 
cesses. Rev. of Geophys. 4 (1), 1-32. 
Hasselmann, K., Barnett, T. P., Bouws, E., Carlson, H., Cartwright, D. E., Enke, 

K., Ewing, J. A., Gienapp, H., Hasselmann, D. E., Kruseman, P., Meerburg, A., 

Muller, P., Olbers, D. J., Richter, K., Sell, W. & Walden, H. 1973 Measurements 

of wind- wave growth and swell decay during the Joint North Sea Wave Project. Deut. 

Hydrogr. Z. 8 (12), 1-95, suppl. A. 
Herbers, T. H. C. & Guza, R. T. 1991 Wind-wave nonlinearity observed at the sea floor. 

part I: forced-wave energy. J. Phys. Oceanogr. 21, 1740-1761. 
Herbers, T. H. C. & Guza, R. T. 1994 Nonlinear wave interactions and high-frequency 

seafloor pressure. J. Geophys. Res. 99, 10035-10048. 
Herbers, T. H. C, Lowe, R. L. & Guza, R. T. 1992 Field observations of orbital velocities 

and pressure in weakly nonlinear surface gravity waves. J. Fluid Mech. 245, 413-435. 
Hillers, G., Graham, N., Campillo, M., Kedar, S., Landes, M. & Shapiro, N. 2012 

Global oceanic microseism sources as seen by seismic arrays and predicted by wave action 

models. Geochemistry Geophysics Geosystems 115, B05302. 
Hughes, B. 1976 Estimates of underwater sound (and infrasound) produced by nonlinearly 

interacting ocean waves. J. Acoust. Soc. Amer. 60 (5), 1032-1039. 
Kanamori, H. & Given, J. W. 1981 Use of long-period surface waves for rapid determination 

of earthquake-source parameters. Phys. Earth and Planet. Int. 27 (1), 8-31. 
Kedar, S., Longuet-Higgins, M., Graham, Frank W. N., Clayton, Robert & Jones, 

Cathleen 2008 The origin of deep ocean microseisms in the north Atlantic ocean. Proc. 

Roy. Soc. Lond. A pp. 1-35. 



double-frequency noise generation by surface gravity waves waves 



33 



Kibblewhite, A. C. & EwANS, K. C. 1985 Wave-wave interactions, microseisms, and infrasonic 
ambient noise in the ocean. J. Acoust. Soc. Amer. 76 (3), 981-994. 

Knudsen, V. O., Alford, R. S. & Emling, J. W. 1948 Underwater ambient noise. J. Mar. 
Res. 7, 410-429. 

Koper, K. D., Seats, K. & Benz, H. 2010 On the composition of Earth's short-period seismic 
noise field. Bull. Seismol. Soc. Am. 100 (2), 606-617. 

Kurrle, D. & Widmer-Schnidrig, R. 2008 The horizontal hum of the Earth: A global back- 
ground of spheroidal and toroidal modes. Geophys. Res. Lett. 35 (2), L06304. 

Lamb, H. 1932 Hydrodynamics, 6th edn. Cambridge, England: Cambridge University Press. 

Landes, M., Hubans, F., Shapiro, N. M., Paul, A. & Campillo, M. 2010 Origin of deep 
ocean microseisms by using teleseismic body waves. J. Geophys. Res. 115, B05302. 

de Laplace, P. S. 1776 Suite des recherches sur plusieurs points du systeme du monde (XXV- 
XXVII). Mem. Presentes Acad. R. Set. Inst. France pp. 542-552. 

Latham, G. V. & Sutton, G. H. 1966 Seismic measurements on the ocean floor. J. Geophys. 
Res. 71 (10), 2545-2573. 

Lighthill, James 1978 Waves in fluids. Cambridge, United Kingdom: Cambridge University 
Press, 504 p. 

Lighthill, J. M. 1952 On sound generated aerodynamically. Proc. Roy. Soc. Lond. A 211, 
564-587. 

Lloyd, Stuart P. 1981 Underwater sound from surface waves according to the Lighthill-Ribner 

theory. J. Acoust. Soc. Amer. 69 (2), 425-435. 
Longuet-Higgins, M. S. 1950 A theory of the origin of microseisms. Phil. Trans. Roy. Soc. 

London A 243, 1-35. 

Longuet-Higgins, M. S. 1970 Mass transport in the boundary layer at a free oscillating surface. 

J. Fluid Mech. 8, 293-306. 
Longuet-Higgins, M. S. & Stewart, R. W. 1962 Radiation stresses and mass transport in 

surface gravity waves with application to 'surf beats'. J. Fluid Mech. 13, 481-504. 
Miller, G. F. & Pursey, H. 1955 On the partition of energy between elastic waves in a 

semi-infinite solid. Proc. Roy. Soc. Lond. A 233 (1192), 55-69. 
Nishida, K. & Fukao, Y. 2006 Source distribution of earth's background free oscillations. J. 

Geophys. Res. 112 (10), B06306. 
Nishida, Kiwamu, Kawakatsu, Hitoshi, Fukao, Yoshio & Obara, Kazushige 2008 Back- 
ground Love and Rayleigh waves simultaneously generated at the Pacific ocean floors. 

Geophys. Res. Lett. 35 (2), L16307. 
Obrebski, Mathias, Ardhuin, Fabrice, Stutzmann, Eleonore & Schimmel, Martin 2012 

How moderate sea states can generate loud seismic noise in the deep ocean. Geophys. Res. 

Lett. 39, L11601. 

Okal, E. A. 1988 Seismic parameters controlling far-field tsunami amplitudes: A review. Natural 
Hazards 1, 67-96. 

Pasyanos, M. E., Walter, W. R. & Matzel, E. M. 2009 A simultaneous multiphase approach 
to determine p-wave and s-wave attenuation of the crust and upper mantle. Bull. Seismol. 
Soc. Am. 99, 3314-3325. 

Peregrine, D. H. 1976 Interaction of water waves and currents. Advances in Applied Mechanics 
16, 9-117. 

Posmentier, E. 1967 A theory of microbaroms. Geophys. J. R. Astron. Soc. 13, 487-501. 
Rhie, J. & Romanowicz, B. 2006 A study of the relation between ocean storms and the Earth-s 

hum. Geochemistry Geophysics Geosysterns 7 (10), Q10004. 
Shapiro, N. M., Campillo, M., Stehly, L. & Ritzwoller, M. H. 2005 High-resolution 

surface-wave tomography from ambient seismic noise. Science 307, 1615-1617. 
Snoke, J. A. 2009 Traveltime tables for iasp91 and akl35. Seismological Research Letters 88 (2), 

260-262. 

Stoneley, R. 1926 The effect of the ocean on Rayleigh waves. Mon. Not. Roy. Astron. Soc. 
Geophys. Suppl. 1, 349-356. 

Tanimoto, T. 2007 Excitation of normal modes by non- linear interaction of ocean waves. Geo- 
phys. J. Int. 168, 571-582. 

Tyler, G. L., Teague, C. C, Stewart, R. H., Peterson, A. M., Munk, W. H. & Joy, 



34 



F. Ardhuin and T. H. C. Berbers 



J. W. 1974 Wave directional spectra from synthetic aperture observations of radio scatter. 

Deep Sea Res. 21, 989-1016. 
Vinnik, L. P. 1973 Sources of microseismic P waves. Pure Appl. Geophys. 103 (1), 282-289. 
Waxler, Roger & Gilbert, Kenneth E. 2006 The radiation of atmospheric microbaroms by 

ocean waves. J. Acoust. Soc. Amer. 119, 2651-2664. 
Webb, S. C. 1998 Broadband seismology and noise under the ocean. Rev. of Geophys. 36, 

105-142. 

Webb, S. C. 2007 The earth's 'hum' is driven by ocean waves over the continental shelves. 
Nature 445, 754-756. 



